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Adiabatic time evolution of degenerate eigenstates of a quantum system provides a means for 
controlling electronic states since mixing between degenerate levels generates a matrix Berry phase. 
In the presence of spin-orbit coupling in n-type semiconductor quantum dots the electron Hamil- 
tonian is invariant under time reversal operation and the many body groundstate may be doubly 
degenerate. This double degeneracy can generate non-Abelian vector potentials when odd number 
of electrons are present. We find that the antisymmetry of many electron wavefunction has no effect 
on the matrix Berry phase. We have derived equations that allow one to investigate the effect of 
electron correlations by expressing the non-Abelian vector potentials for many electron system in 
terms of single electron non-Abelian vector potentials. 
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I. INTRODUCTION 



Single electron control in semiconductor quantum dots 
would be valuable for spintronics, quantum informa- 
tion, and spin qubits[l], Adiabatic time evolution 
of degenerate eigenstates of a quantum system provides 
a means for controlling individual electrons since mix- 
ing between degenerate levels generates a matrix Berry 
phase 0, It has been shown that universal quantum 
computation is possible by means of non-Abelian uni- 
tary operations [Bj, [f|. Recently several proposals on how 
to generate matrix Berryphases have been made in semi- 
conductor systems 0, H, M, OH EH- ^ n one °f these pro- 
posals it has been suggested that a single electron in a 
n-type semiconductor quantum dot with spin-orbit terms 
can have matrix Berry phasesfll]. All the discrete energy 
levels of quantum dots possess a double degeneracy be- 
cause the Rashba and/or Dresselhaus spin orbit coupling 
terms have time reversal symmetry jl2j|. In Ref.[ll| it was 
suggested that the adiabatic transformation can be per- 
formed electrically by changing the confinement poten- 
tials of the quantum dot. It was shown that matrix Berry 
phases can be produced only when the parity symmetry 
of two-dimensional harmonic potentials is broken. A de- 
generate pair of states also exists in quantum rings (T^. 
The relevant symmetry of this system is the combined 
operation of time reversal and large gauge transforma- 
tions. However, it is often not easy to populate a quan- 
tum dot or ring with just one electron and usually there 
are several electrons in it [3 EH- It is not certain that 
many body groundstates are doubly degenerate in the 
presence of time reversal symmetry. Moreover, it is un- 
clear how the effects of exchange and correlation many 
body physics change the non-Abelian vector potentials. 
These issues are investigated in the present paper. 

Before we give a summary of the main points of the 



present paper let us give a brief introduction of the basic 
ideas behind matrix Berry phases. Consider a degener- 
ate pair of states and l^), which may be single 
particle or many body states. If the system is in a su- 
perposition state |*(0)) = ci(0)|*i(0)) + c 2 (0)|* 2 (0)) 
at time t=0 an adiabatic evolution of the parameters 
X p can transform this state into another state = 
ci(£)|*i(t)) + c 2 (t)\^ 2 (t)} after some time t. Here the 
orthonormal basis states |$i(t)) are the instantaneous 
eigenstates of the Hamiltonian. For a cyclic change with 
the period T, represented by a closed contour C in the 
parameter space, the states |^ , i(T)) and \^f 2 (T)) return 
to the initial states |*i(0)) and |*2(0)), but the coeffi- 
cients c\(T) and c 2 (T) may not return the initial values. 
In such a case a 2 x 2 matrix Berry phase (non-Abelian 
Berry phase) &c is generated 
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The expansion coefficients c\(t) and c 2 (t) satisfy the 
time-dependent Schrodinger equation 



ihti 



1,2. 



(2) 



The 



matrix elements 



Aij are given by A 



ill 



where the sum over p in Ai 



is meant 



to be the sum over A p . The non-Abelian vector poten- 



tials A p are 2x2 matrices 



(A p )i y: 
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Note that the wavefunctions and tyj are degenerate. 
The role of the non-Abelian vector potentials is to keep 
the energy degeneracy intact while the system parame- 
ters change. 

Our investigation shows that when even number of 
electrons are present the many body groundstate may 
not be doubly degenerate. However, when odd number 
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FIG. 1: Single electron eigenstates form degenerate energy 
shells. In each shell electron wavefunctions form a time re- 
versed pair. 



definition of single electron non-Abelian gauge potentials 
is extended to 



(5) 



when (f>k and <f> p belong to different single electron energy 
shells. We will call these inter shell non-Abelian gauge 
potentials. 
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of electrons are present the groundstate may be dou- 
bly degenerate. In the simplest possible approximation 
the antisymmetry of many body wavefunctions can be 
ignored. Then the degenerate groundstates for three 
electrons are given by \$>i) = (j)i(ri)(/)2(r2)(j>3(r3) and 
|<J7i) = 0i(ri)02(r2)</ , 4(^3), where 4>i are single electron 
wavefunctions, see Fig[T] In this case one can show that 
the solutions of the time dependent Schroedinger equa- 
tion, Eq.([2]), have the form 

\c 2 {t)) \B 2 {t))- ^ 

The function f{f) depends on the single electron non- 
Abelian vector potentials of the first shell, and the func- 
tions B\ (t) and B 2 (t) depend on the single electron non- 
Abelian vector potentials of the second shell, see Figfj] 
(The function f(t) is periodic with the period T ). The 
probability to find an electron in the state | <^>3 ) ( 1 04 ) ) is 
| Cl (T)| 2 = |i? 1 (T)| 2 (|c 2 (T)| 2 = |£? 2 (T)| 2 ). So these prob- 
abilities are independent of the phase factor e^A*), i.e., 
of the non-Abelian vector potentials of the first shell. 

In the next simplest possible approximation the anti- 
symmetry of many body wavefunctions can be included. 
In this case the doubly degenerate many body ground- 
states may be taken to be two Slater determinant wave- 
functions, and ^i, that are time reversed states of 
each other, see FigEJ It will be shown below that the 
antisymmetry of these Slater determinant wavefunctions 
does not change the results obtained above for the ma- 
trix Berry phase. In addition to the antisymmetry one 
can include correlation effects. This can be done by in- 
cluding four basis vectors \&i ^2 and ^2, see Figj2] 
In this case the groundstates are given by linear combi- 
nations of these four states and are thus correlated. It 
can be shown that the non-Abelian gauge potentials for 
a correlated electron system are related to non-Abelian 
gauge potentials for single electrons, provided that the 



FIG. 2: Four many body basis states in the three electron 
Hilbert space. The basis pair ^Pi) or (^J/i , ) is a time 
reversed pair and has the same total single electron energy. 
So do V&2 and ^2- Black dots represent electrons. 



In Sec. II we describe the Hamiltonian of the sys- 
tem in detail, and in Sec. Ill we discuss why matrix 
Berry phases can be produced only when the parity sym- 
metry of two-dimensional harmonic potentials is broken. 
A computational scheme to include systematically many 
electron physics is explained in Sec. IV. The absence of a 
matrix Berry phase for the two-electron case is shown in 
Sec. V. In Sec. VI we show that a matrix Berry phase is 
present in the three-electron case. Correlation effects are 
discussed in VII and conclusions are given in Sec. VIII. 

II. MODEL HAMILTONIAN 

An electric field E is along the z-axis and electrons 
are confined in a triangular potential V(z). When the 
the width of the quantum well along the z-axis is suffi- 
ciently small we may include only the lowest subband 
state along the z-axis. We denote this wavefunction 
by f{z). From the expectation value (f{z)\k 2 \f{z)) = 
0.8(27n* eE I ft 2 ) 2 / 3 we estimate the characteristic length 
scale along the z-axis: R z = 1/ ^0.8(2m*eE/h 2 ) 2 / 3 . The 
Hamiltonian in the absence of the spin orbit coupling is 

H K = -^ + U(x,y) + V(z). (6) 

We take the two-dimensional potential to be U(x,y) = 
^m*Lu 2 x 2 + ^m*Lu 2 y 2 +V p (x, y). The strengths of the har- 
monic potentials are denoted by uj x and u) y . They may 
be varied by changing gate potentials of the quantum dot 
system. The characteristic lengths scales along x- and 
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y-axis are R x v — \/t^7 l , — ■ The potential V v (x,y) is 

chosen in such a way that the parity symmetry of elec- 
tron wavefunctions is broken. We can choose it to be 
V p (x,y) — e'y. It represents an electric field along the 
y-axis and its strength e' may be varied electrically. 

In a periodic crystal potential of a semiconductor the 
spin orbit interaction has two contributions. The Rashba 
spin orbit termfl6j is 



The single electron wavefunctions used in this section 
are self-consistent Hartree-Fock single electron states if 
Hs is replaced by the HF single electron Hamiltonian. 
In the absence of the Zeeman term the single particle 
Hamiltonian is invariant under time reversal symmetry: 
k — > —k and & — ► — a. The time reversal operator is 
T = —icryC, where the operator C stands for complex 
conjugation. The time reversed state of \<j>) is 
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(7yk x 



(7) 
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\mn t) +^< mT |mn J.) {11) 



Here a x , y are Pauli spin matrices and k x , v are momentum This can also be written as a co l um n vector |&) 



operators (k x = and similarly with k y .). The con- 
stant c R depends on the external electric field E applied 
along the z-axis. The Dresselhaus spin orbit term[17| is 

H B = cd {(a x k x (k 2 y - k 2 z )) + (a y k v (k 2 z - k 2 x ))) . (8) 

There is another term of the form a z (k z ) (k 2 . — k 2 ) in 
the Dresselhaus spin orbit term but it vanishes since 
the expectation value (k z ) = (f(z)\k z \f(z)) = for the 
first subband along z-axis. The constant cjj represents 
breaking of inversion symmetry by the crystal in zinc 
blende structures. The total single electron Hamiltonian 
of semiconductor quantum dot is Hs = i?K + -Hr + Hd ■ 
When several electrons are present in the dot the to- 
tal Hamiltonian consists of single particle Hamiltoni- 
ans and two-particle interaction terms H = H + Vi nt , 
where Hq 

V(n 



J2t H t and V int = \ Y,i^j V {^i ~ *j) witn 
{Hi is the single particle Hamilto- 
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nian for the ith electron and e is the background dielec- 
tric constant). The total many electron Hamiltonian is 
invariant under time reversal symmetry. 



III. BREAKING OF PARITY SYMMETRY OF 
WAVEFUNCTIONS 

The eigenenergies of the single electron Hamiltonian is 
depicted in FigJTJ Each energy is doubly degenerate. An 
eigenstate of the Hamiltonian Hs may be expanded as 
a linear combination of the eigenstates of Hk ■ A lowest 
energy state can be written as 



|0i 



y^cwrlmn T) + y^c OTW ;lwn I)- (9) 



In the basis states \mno~) the quantum number m(n) and 
a label the harmonic oscillator levels along the x-axis (y- 
axis) and the component of electron spin. The subband 
wavefunction f(z) is suppressed in the notation \mna). 
The wavefunction \<f>i) can also be written as a column 

5), where F T = E^-tI™) and 

F[ = J2mn c mni\rrm). The expansion coefficients satisfy 
a matrix equation 



vector \<p\ 



(mna\H\m'n'a') 



= Ec„ 



(10) 



Ft 



Note that T 2 | 



These two states 



are degenerate and orthonormal. We have suppressed 
the Bloch wavefunction of the conduction band in ap- 
plying the time reversal operator since it is unaffected 
by the operator T. Our wavefunctions are all effective 
mass wavefunctions and only the conduction band Bloch 
wavefunction at k = is relevant. 

Even in the presence of doubly degenerate states a ma- 
trix Berry phase may not always exist. It was shown nu- 
merically that only distorted two-dimensional harmonic 
potentials breaking the parity symmetry generate matrix 
Berry phases[ll|. Here reasons for this effect for the sin- 
gle electron case is given in detail. In this argument only 
one electron is assumed to be present in the dot, crj = 0, 
and Vp(x, y) = 0. The adiabatic parameters are Ai = cr 
and A2 = u) x . Let us consider Hr as a perturbation and 
apply a perturbation theory to 1 00 f) 



|0> a |00 T) 



\mn i){mn[ \H R \QQ f) 



III II 



EE 



Eqq — E mn 

mn t)(mn T \H R \m'n' |) )<mW { \H R \QQ T) 



no 



E„ 



Eqo — E„ 



+ (12) 

where E nm = (m + l/2)huj x + (n + l/2)hu> v . The wave- 
function of 1 00 t) has even parity. The matrix element 
(mn I \H R \00 f) couples |00 f) only to odd parity states 
\mn I). Therefore the wavefunction corresponding to 
the second term in Eq. (|12[) has odd parity with spin 
down. The wavefunction of the third term in Eq. (fT2"]) 
has even parity with spin up. This perturbation series 
indicates that the exact lowest energy state has the form 

\4>i) = ( _f 1 , where F e and F a even and odd functions. 



F 



The time reversed wavefunction is |^) 



_ p* 

x o 

f; 



For 



this particular pair of degenerate states the off-diagonal 
elements of non-Abelian vector potentials are zero: 



1^- 

9A, 



i,p k q,p k 







(13) 
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Here F e = J2 q a q l Pq with q = (m, n) such that the two- 
dimensional harmonic wavefunctions ip q are even func- 
tions. The function F Q = ^2b p (p p with p = (m',n') are 
such that ip p are odd functions (Note that one can show 
(</3 g | gf-l^p) = 0). Since the off-diagonal non-Abelian 
vector potentials are zero the matrix Berry phase is ab- 
sent. In what follows a distorted two-dimensional har- 
monic potential is assumed. The adiabatic parameters 
are cr and e'. 



IV. COMPUTATIONAL SCHEME TO INCLUDE 
MANY ELECTRON PHYSICS 



E 



¥1 ¥ 2 ¥ 2 ¥ 3 ¥, 

FIG. 3: Five many body basis states in the two electron 
Hilbert space. The basis state ^1 is the time reversed state 
of itself. 

A computational scheme that includes systematically 
many electron physics is explained here. We find an ef- 
ficient computational scheme of computing eigenstates 
and eigenvalues of the many electron Hamiltonian of a 
semiconductor quantum dots with an arbitrary shape of 
the confinement potential. We follow the following steps. 

Single electron Hamiltonian is first solved in the pres- 
ence of spin-orbit terms. This can be done by following 
the steps given in Sec. III. The number of single electron 
basis states \mna) can be truncated to a reasonable num- 
ber by comparing the relative magnitude of the single 
particle energy spacing and the strength of the Coulomb 
interaction. 

Using these single electron eigenstates Slater deter- 
minant states and their time reserved states are 
formed. ^ is obtained from by replacing the single 
electron wavefunctions with their time reversed states: 
<pk — > <p k . It will be inefficient to use two dimen- 
sional harmonic potential wavefunction ip mn to build 
these Slater determinant states. The following relation 
hold for two arbitrary states ^ a and * b : (T^ a \T^ b ) = 
(^bl^a). From this relation it follows that the Hamil- 
tonian matrix elements have the following properties: 
{*i\H\* t ) = (9i\H\*i), (*i\H\*j) = -i^fjlH^i), and 
(%\H\^) = (*i\H\9j)*. 

When there are odd number of electrons the basis 
vectors of many electron Hamiltonian are chosen to in- 
clude time reversed pairs: {"Fi, <Fi, ^ 2 , ^2, ■•■}, see Fig. 
[3J The eigenstates are doubly degenerate in the odd 
electron case. Let us show this for the three electron 



case using the following Slater determinant wavefunc- 
tions as the basis vectors in the many body Hilbert 
space {^1,^1,^2)^2}' These states are depicted in 
Figl^l The wavefunctions of the lowest energy degener- 
ate pair are $ = a*i + 6*1 + c* 2 + d^ 2 and ¥ = T$ = 
a*Wi-6**i + c** 2 -d**2. Clearly {$\H\®) = 
Note $ and $ are a time reversed pair, orthogonal, and 
satisfies T 2 <E> = — $. The two eigenvalues of H are each 
doubly degenerate. 

When there are even number of electrons the basis 
vectors of the many electron Hamiltonian are chosen as 
{^1, ^2, ^2, ^3, ^3, •■•}• In this scheme the groundstate 
is non degenerate. This can be shown shown, for exam- 
ple, in an approximation that includes five many body 
basis vectors in the presence of two electrons, see Fig. 
121 The groundstate can then be written as a linear com- 
bination $ = a*i + 6*2 + c¥ 2 + d^ 3 + e* 3 . Its time 
reversed state is ¥ = a**! + b*H 2 - c** 2 + d**3 - e*^ 3 . 
Clearly these two states are not orthogonal. Since 
= (T 2 $\H\T 2 $) and T 2 |<F) ^ -|$) one finds 
that the energy expectation values are non degenerate, 
^ i.e., time reversed states are non 

degenerate. The matrix Berry phase is thus absent. In 
the simplest possible approximation the groundstate 'F 
of two electron system is given by the single Slater de- 
terminant made of the single particles states {0i,</>i}. In 
this case the initial state should return to itself after the 
adiabatic cycle. It is instructive to verify this explicitly 
by transforming each single particles states. The initial 
state is 

|#(0)) = <Mri)^(r 2 ) - Mr2)Mn). (14) 

After a cycle the single electron wavefunctions transform 
into 

4>i — > ai</>i + b\<j) 1 and <j) 1 — ► a 2 (j)i + 62^1 ■ (15) 

The transformed two electron Slater determinant is 
|*(T)) = (016a - ha 2 )(Mn)Mr2) - 0i(Wi(ft)), 
which is (<ii&2 — &i«2)|^(0))- Since the transformed sin- 
gle electron states are time reversed states of each other 
(02,62) = ( — b*,a*). So the factor in the transformed 
wavefunction is aib 2 — b\a 2 = \ci\\ 2 + |6i| 2 = 1. 



V. SIMPLEST POSSIBLE APPROXIMATION 
WITH ANTISYMMETRY OF WAVEFUNCTIONS 

Now consider odd number of electrons, for example, 
three electrons in a dot. In the simplest possible approx- 
imation with fermionic antisymmetry two many body 
basis states, <Fi and <3>i, can be included. These are 
Slater determinant states and are depicted in Fig[2j This 
approximation is well justified when the typical single 
electron energy spacing is larger than the characteris- 
tic Coulomb energy. Clearly 'Fi and VPi are orthogonal 



5 



(^il^i) = 0. The Hamiltonian matrix is 
H = 



(£i|ff|*i) 



£i|ff|£i) 
*i|ff|*i> 



The coefficients c\(t) and C2(i) satisfy the Schroedinger 
equation Eq.([2]) governed by the non-Abelian vector po- 
tentials 



(16) 



The off-diagonal matrix elements of the Hamiltonian ma- 
trix is zero since {^i\H\^i) = -(<I' 1 |i7|l' 1 ), which fol- 
lows from (T* a |T* 6 ) = (* 6 |* a ). It is instructive to 
show this explicitly by evaluating the matrix element 

= [(23|F|24) - (23|V|42)] + [(13|F|14) - (13|V|41)]. 

(17) 

Here the first and third terms are direct Coulomb ma- 
trix elements while the second and fourth terms are their 
exchange counter parts. The direct Coulomb matrix el- 
ements are zero since 4> 3 (r)4>^(f) = (j) 3 (r)(f> 3 (f) = 0. So 
only their exchange counter parts remain, which are given 
by 



(^fe)2,2 



(24) 



/ 



It can be shown using Eqs. (|A3p and (|A4[) that these can 
be written in terms of the non Abelian vector potentials 
for single particle states 

(Ak)l,l = (afc)l,l + (afc)2,2 + (Ct/s)3,3) 

(Ak)l,2 = (afc)3,4, 

(Ak)2,l — (afe)4,3, 

(A k )2,2 = (afe)i,i + (afe)2,2 + (afc)4,4, (25) 



(23|y|42> = / dndriMnWs^Win - r 2 )0 3 (r 1 )0 1 (r 2 ) 



(18) 



where the non-Abelian vector potentials involving the 
lowest energy single electron states are 



and 



(13|V|41) = / dr 1 dr 2 <t>t(?i)cf>$(r 2 )V(r 1 -f2)<f> 3 (n)Mr f 2). 



(Ofc)l,l : 

(0^)2,2 



(19) 



O* 



(26) 



and involving the second lowest single electron states are 



One can show 

(^i\V int \^i) = -<23|y|42) - <13|V|41> = 0. 
This is because 

<23|V|42) = — <13|V|41>, 



(20) 



(21) 



which follows from 

0* (r) 03(f) = 
^(^(r) = 

0!(r> 3 (r) = 



J 3- 



and 



(ctfe) 


'3,3 = *( 




(flfe) 


U,4 = 






! 3,4 = i( 




(Ofc 


U,3 = 





(27) 



The solution of the time dependent Schroedinger equa- 



tion has the form 



B x (t) 



The func- 



(22) tion A(i) depends on the electrons in the first shell. It 



Note that the second and third terms are equal except 
for the minus sign. The off-diagonal elements are thus 
zero, and 



satisfies f(t)A(t) = -ifh%, where f(t) = Efe(K)i,i 
(0/0)2,2)^^- The solution of this differential equation is 
a phase factor 



H 











A(t) 



.*/(*) 



(28) 



(23) 



This implies that \&i and 'J'i are eigenstates and are de- 
generate in energy. 

During the adiabatic change the wavefunction of three 
electrons evolves as |*'(i)) = ci(*)|*i(*)) + c 2 (t)|*i(t)). 



Note that (a,k)i,i + (0^)2,2 appears in the diagonal ele- 
ments of the non Abelian vector potentials (Ak)n and 
(Ak)2,2- The functions B±(t) and B2(t) are the solutions 
of the time dependent Schroedinger equation 



ihBi 



1,2. 



(29) 



with the following non-Abelian gauge potentials. 



The non-Abelian vector potentials are 





= (Ofc)3,3) 




(Ak)l,2 


= (Ofc)3,4, 






= (afc)4,3j 




(A k ) 2 ,2 


= (0/0)4,4- 


(30) 



This is precisely the non Abelian vector potentials for the 
second energy shell. The probability to find an electron 
in the state |<fe)(|<M) is \ci(T)\ 2 = |Bx(T)| 2 (|c 2 (T)| 2 = 
|-B 2 (T)| 2 ). So these probabilities are independent of the 
non-Abelian vector potentials of the first energy shell. In 
addition it should be noted that the antisymmetry of the 
many electron wavefunction does not change the value of 
the matrix Berry phase obtained from the wavefunctions 
without the antisymmetry, see Eq.(|3]). This is the main 
conclusion of this section. 



VI. CORRELATION EFFECTS 

Correlation effects can be taken into account by writing 
many body eigenstates as a linear combination of single 
Slater determinant wavefunctions {<E>i, <I> 2 , $3, 3>4, ■■■} = 
{^1, ^1, ^2, ^2, ■■■}• These wavefunctions arranged in 
the increasing order of the total kinetic energy E\ < 
E 2 < .... Note that ($i|<l?j) = 5 t j for all i and j. If the 
groundstate is doubly degenerate one of the degenerate 
groundstates can be written as 



A I 



|$)=X)ci|#i), 



(31) 



with 



A I 



Y,{HE\*j)ci = Eci 
3=1 



(32) 



Then its time reversed state has the form 

At 



|$>=T$ = ^ c *|$ l 



(33) 



i=l 



since T is an antiunitary operator. Note that 
<&j are ordered as follows {$1, $ 2 , <j> 3 , $4, ...} = 
{Wi,*xi*2)*2, — }■ The states |4>) and |¥) are eigen- 
states of the truncated Hamiltonian with the dimension 
M. (Note that as M — » 00 the energy E will approach 
the exact eigenvalue.). Some of the Slater determinant 
states are depicted in Fig[2j Clearly |$) and |$) are 
orthogonal. The matrix element between them is zero: 
($|iJ|¥) = 2?($|$)=0. These two states are degenerate 
since {$\H\$) = @\H\W). 



(A k )x,i 


















( A kh,2 




(34) 



In terms of many body basis states they can be written 
as 



(A k )i,i 
(Akh,2 

(Ak)l,2 

(Akh,i 



Q c * Q 



d ,_ 



(35) 



In the three electron case with the basis vectors given 
in FigJ5] the diagonal matrix elements z($j|g^-|<l>j) and 

*( $ il5irl $ ») are S iven b y 



(flfc)i,i + (0*3)2,2 + (ofc)3,3j 
+ (afc)a,2 + (afc)4,4, 
(afe)i,i + (ofc)a,2 + (0*1)5,5, 



i<*x|^|*x) 
^1^1**} 

«^l*2) 

— d — 

*<*2 kH^a) = (a fc )i,i + (a k ) 2 ,2 + (afc)6,6- (36) 

Similarly one can find off-diagonal matrix elements 
d 



_9_ 

9A fc 

9A fc 



l*i> 
1*2) 
1*2) 
|W 2 ) 
l*2> 
1*2) 



( a fe)3,4, 

(a*)s,6) 
(0*1)3,5, 
(0*1)3,6, 

(Ofe)4,5, 
(ofc)4,6- 



(37) 



These results follow from Eq. (IA4j) . Usually the non 
Abelian vector potentials are defined between two de- 
generate states, see Eq.Q. When electron correlations 
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are important the matrix elements between two single 
electron states belonging to different energy shells are 
relevant (see the last four equations in Eq. ([37|) ). We will 
call these inter shell non-Abelian vector potentials. Phys- 
ically the mixing between the second and third shells is 
facilitated by them. Note that {ak)j,% = (a>k)*j- 

VII. DISCUSSIONS AND CONCLUSIONS 

We have developed a scheme of computing eigenstates 
and eigenvalues of the many electron Hamiltonian of a 
semiconductor quantum dot with an arbitrary shape of 
the confinement potential. The basis vectors of this many 
electron Hamiltonian are chosen to include time reversed 
pairs of Slater determinant states. In this scheme the 
groundstate is non degenerate when an even number of 
electrons are present while the eigenstates are doubly de- 
generate when an odd number of electrons are present. 

Using this computational scheme we have investigated 
matrix Berry phases when several electrons are present 
in a semiconductor quantum dot with spin-orbit coupling 
terms. When an odd number of electrons is present a ma- 
trix Berry phase exists. The effect of the antisymmetry 
of many electron wavefunctions does not affect the value 
of the matrix Berry phase. Electron correlation effects 
can be treated by including the inter shell coupling. In 
this case the many body non-Abelian vector potentials 
can be written as a sum of those of the single electron 
non-Abelian vector potentials, provided that the inter 
shell non-Abelian vector potentials, i{4>p\-§\\(t>q) with <j> p 
and 4> q non degenerate, are introduced to represent the 
shell coupling. It would be interesting to investigate the 
role of these inter shell non-Abelian vector potentials in 
parabolic dots. Optical dipole transitions (l8l. [l9| can be 
used to detect matrix Berry phases even when several 
electrons interact strongly with each other. This is be- 
cause the dipole oscillator strength will depend on the 
matrix Berry phasefll]. However, our results in Sec. V 
indicate the matrix Berry phases will be independent of 
many body physics when the energy level spacing be- 
tween the shells is larger than the characteristic Coulomb 
energy. In self assembled quantum dots with wetting lay- 
ers this condition is satisfied [20|. In this case only the in- 
tra shell non-Abelian vector potentials, i{(t>p\-§x\4>q) with 
(j> p and 4>q degenerate, are relevant, and the value of the 
matrix Berry phase can be calculated as in ref.[ll[. We 
leave it as a future study to investigate numerical val- 
ues of the matrix Berry phases for strongly interacting 
electrons. 

In II- VI semiconductors the Rashba term is expected 
to be larger than the the Dresselhaus coupling and in 



III-V semiconductors the opposite is true[16|. The typi- 
cal value of the energy scale associated with the Rashba 
constant depends on the electric field applied along the 
z-axis and the semiconductor material: it is of order 
Er = cr/R ~ 0.01 — lOmeV, where the length scale 
R ~ lOOA is the lateral dimension of the quantum 
dot. The second adiabatic parameter E p represents the 
strength of the distortion potential e'y: the expectation 
value of the distortion potential is E p — (0\e'y\l). Its 
magnitude is of order 1 — lOmeV, depending on the elec- 
tric field applied along the y-axis. 

In self-assembled InAs quantum dots [20] the level spac- 
ing can be larger than the characteristic Coulomb en- 
ergy, and our investigation shows that in such a case 
the many problem essentially reduces to a single parti- 
cle problem (The energy E associated with the length 
scale R is typically of order 20-40meV in self-assembled 
quantum dots). In this regime, as we explained above, 
the exchange effect disappears in matrix Berry phases 
and it is only the last singly occupied electron that mat- 
ters. The presence of the other electrons in the inert 
filled shells will just modify the confinement potential. 
Thus the matrix Berry phase can be calculated follow- 
ing the method outlined in Ref[ll], and below we give 
some typical values of the matrix Berry phases. Wo 
take an elliptic adiabatic path given by (Eii(t), E p (t)) = 
{Er, c + AEr cos (wt), E PtC + AE P sin(wt)). For the pa- 
rameters Er^/Eq = 2, E p , c /E = 1, AE p /Eo = 
0.9, and hLj/E = 0.1 we find (\ Cl (T)\ 2 ,\c 2 {T)\ 2 ) = 
(0.7971,0.2036), (0.9138,0.0867), (0.9810,0.0191), re- 
spectively, for different values of the semiaxis of the 
ellipse AEr/Eq — 1.9, 1.5, 1 (The initial state is 
(|ci(0)| 2 , |c 2 (0)| 2 ) = (1,0)). Note that the parameter for 
the semiaxis AEr/Eq can be varied experimentally ei- 
ther by changing the Rashba field AEr or the size of 
the dot through Eq. We have verified that two differ- 
ent values of the frequency uj give the same values of 
(|ci(T)| 2 , |c2(T)| 2 ) as long as the other parameters are 
unchanged, confirming that the matrix Berry phase is 
a geometric phase [3] ■ This should be tested in any ex- 
perimental investigation of matrix Berry phases. Exper- 
imental data of single quantum dots would be less com- 
plicated than data of an ensemble of dots with various 
sizes of dots, and can thus be compared directly with our 
theoretical results. 



APPENDIX A: MATHEMATICAL NOTE 

Let consists of {<f>i, <j>2, 03 } and let \^2) consists 
of {0i, 02, ^4}: 
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l*i> = 



V3! 



{ 0i(ri)02(r 2 )0 3 (r 3 ) + 03(ri)0i(r2)02(r 3 ) + 02(ri)03(r2)0i(r 3 ) 
- 03(ri)02(^ 2 )0i(r3) - 0i(ri)03(r 2 )02(? , 3) - 02(ri)0i(r 2 )03(r-3)} 



(Al) 



and 



|*2> = -j= { 0i(ri)02(r2)04(7-3) + 04(ri)0i(r2)02(r3) + 02(ri)04(r 2 )0i(r 3 ) 

- 04(ri)02(r 2 )0i(r 3 ) - 0i(ri)04(r 2 )02(r 3 ) - 02(ri)0i(r 2 )04(r 3 )}. 

I 



(A2) 



The diagonal elements of the non-Abelian gauge poten- 
tial between the Slater determinant wavcfunctions are 



Y.'MrH^-JM''))- (A3) 



If more than two single electron wavefunctions in 
must be exchanged with those in to transform 
into |*,) then (^| gf-|*,) = 0. 



Let us find the off-diagonal non-Abelian vector potential 
elements i(^ r i|g^-|^ r 2) between the two Slater determi- 
nant wavefunctions |^i) and \^f 2 ) ■ When 3 is replaced 
with 04 the state |<3/i) transforms into \^2)- We find 
(*ilsjbl*2> = (03(^)1^104^)). This result may be 
easily generalized to N particle Slater determinant wave- 
functions. In general the non-Abelian vector potential 
element between two different Slater determinant wave- 
functions, \^fi) and is non zero when can be 
transformed into by replacing one single electron 
wavefunction in l^) with another wavefunction in 
If these single electron wavefunctions are denoted by 0£ 
with 0„ then we find 



mi-g^iM?))- (A4) 
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